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Abstract 

We consider N-body problems with homogeneous potential 1/r^” where 
K £ (0, 1), including the Newtonian case (k = 1/2). Given R > 0 and 
T > 0, we find a uniform upper bound for the minimal action of paths 
binding in time T any two configurations which are contained in some 
ball of radius R. Using cluster partitions, we obtain from these estimates 
Holder regularity of the critical action potential (i.e. of the minimal action 
of paths binding in free time two configurations). As an application, we 
establish the weak KAM theorem for these N-body problems, i.e. we prove 
the existence of fixed points of the Lax-Oleinik semigroup and we show 
that they are global viscosity solutions of the corresponding Hamilton- 
Jacobi equation. We also prove that there are invariant solutions for the 
action of isometries on the configuration space. 


1 Introduction 

Let E he a, finite dimensional Euclidian space, and denote by a; = (ri,..., tat) S 
the configuration vector of N punctual masses mi,..., uin > 0. By || x || we 
will denote the norm given by max { II n\\E I 1 < i < N }, and | x \ will denote 
the norm induced by the mass scalar product 

N 

<x,y >= ^ mj < Tj, Si >E 

i=l 

for X = (n,..., Tat), y = (si,...,SAr) G E^. As usual, we call I{x) = | a; p 
the moment of inertia of x regarding the origin of E. The A-body problem is 
determined once the force function U on E^ (or potential function), negative 
of the potential energy, is chosen. In this paper, we restrict us to the potential 
functions which are homogeneous of degree —2k 

Uk{x) = ^ miTTlj {Vij) , 
i<j 

where = || — rj || £;, and k G (0,1). The case k = 1/2 corresponds to 

the Newtonian potential. In other words, this means that the laws of motion 
are given on the open and dense subset fl = { a: £ E^ \ Uk.{x) < -boo } by the 
differential equation x = Vf7„, where the gradient is taken with respect to the 
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mass scalar product on . The equivalent variational formulation is given by 
the Lagrangian defined on TE^ = E^ x E^, 

1 ^ 

L{x, v) = rriivf + U^{x), 

^ i=l 

where v = (vi,..., u„). Thus, motions are characterized as critical points of the 
Lagrangian action ^(7) = f ^(7(5), 7(5)) ds, and the Euler-Lagrange equations 
define a - non complete - analytical flow on the non compact manifold Tfl. 

1.1 Globally minimizing curves and the action potential. 

Let us give a precise definition of the Lagrangian action functional. Recall that 
a curve 7 : [a,b] E^ is absolutely continuous if it is differentiable almost 
everywhere, and its derivative 7 satisfies the fundamental theorem of calculus 
for the Lebesgue integral. Thus the Lagrangian action is well defined on the 
set of absolutely continuous curves C. More precisely, the action is the function 
A : C ^ ( 0 , +00] given by 

^(7) = y L{-f{s),'j{s))ds = ^ J \'j{s)\'^ds + J C/«,(7(s)) ds . 

where | z; | is the norm in E^ induced by the mass scalar product. It can be seen 
that absolutely continuous curves with finite action are necessarily l/ 2 -H 61 der 
continuous, hence they are contained in the Sobolev space H^{[a,b],E^). 

For T > 0 and x,y G E^, denote by C{x,y,T) the set of all absolutely 
continuous curves 7 : [ 0 ,T] E^ which satisfy 7(0) = x and 7(T) = y. We 
are interested in the function (p defined on E^ x E^ x (0, +00) by 

y, T) = inf { ^(7) | 7 G C(a:, j/, T) } . 

We will say that a curve 7 : [a, b] —>■ E^ is globally minimizing, if we have 
that 71(7) = <p{pi{a), 7(6), b — a). For a curve defined on a non compact interval, 
globally minimizing will mean that the property is satisfied for all restrictions 
of the curve to a compact interval. It is not difficult to see that a globally 
minimizing curve always exists for any two configurations x,y € E^ and for all 
T > 0 . Essentially, it is a consequence of the lower semi-continuity of the action 
functional. 

In the last years, the global variational methods have been successful to 
prove the existence of a great variety of particular motions. A typical example 
is the eight choreography of Chenciner and Montgomery [ 3 ] , among many others 
closed orbits with topological or symmetry constraints. The main difficulty that 
raises from these methods for the Newtonian potential, and also for the homo¬ 
geneous potentials here considered, is the one to assure that global minimizers 
avoid collisions, that is to say, that they are contained in the open domain 17 . 
Following an idea of Marchal, Chenciner established a proof of this fact, for 
the Newtonian N-body problem in the plane or the three-dimensional space, 
see 0 ) [IS]- Simultaneously and independently, Ferrario and Terracini gave an 
improved version of the Marchal’s theorem, see [ 12 ]. We will nowhere use this 
result in this paper, but it is relevant to remark that combined with proposition 
dH below, and using results from | 7 ], we can deduce (for the Newtonian case in 
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dimension greater than one) the existence of completely parabolic motions with 
arbitrary initial configuration. Recently, this last result was improved in M- 
Our first result gives an upper bound for the action of such curves which 
depends on the size of the configurations. In our opinion, this result is quite 
fundamental for global variational methods, and it is optimal, in the sense that 
the bound is reached by homothetic minimizing configurations, as we explain in 
the following section. 

Theorem 1. There are positive constants a, /3 > 0 such that for all T > 0, 
(fix,y,T) < aT-^R^ + PTR-'^^ , 

whenever x and y are configurations contained in a ball of radius R > 0 of E. 
The constants a and (3 only depend on the degree of homogeneity of the potential 
(—2k), the number of bodies N, and their masses. 

The next result shall be useful for the study of free time minimizers, that 
is to say, absolutely continuous curves which minimizes the action in the set of 
curves C{x,y) = [Jj,^QC{x,y,T). The Mane’s critical action potential (see for 
instance 0), or the action potential, is defined in our setting on x E^ by 

fiix, y) = inf { (j){x, y,T) \ T > 0 } = inf { A{-f) \ 7 e C(a;, j/) } . 

It is clear that 4>(x,y) = (j){y,x), and that 4){x,y) < (j){x,z) + (j){z,y), for any 
conhgurations x, y, z in E^. In fact, proposition [5] shows that the action po¬ 
tential (j> is a distance function. Notice that as a corollary of theorem [U we have 
that (j){x, y) < {a + I3)R^~^ whenever x and y are configurations contained in a 
ball of radius R > 0 of E. With similar arguments as in theorem [TJ combined 
with a cluster decomposition, we obtain the following theorem. 

Theorem 2. There is a positive constant rj > 0 such that for all x,y € E^, 

4‘{x,y) <r]\\x-y\\^~'^. 

Therefore, the action potential is Holder continuous respect to the Euclidean 
norm on E^ X E^. In other words, for any configurations x, y, z in E^ we 
have (j){x, z) — (j){y, z) < (j){x, y) < r]\\x — y\\ On the other hand, it is easy to 
prove that the action potential is locally Lipschitz in the open and dense subset 

nxflC E^ X E^. 

The action potential (j) was introduced by Mane in the nineties, as well 
as the Fathi’s weak KAM theorem, for the study of the dynamics of Tonelli 
Lagrangians on compact manifolds. But in fact, for the Newtonian case of our 
N-body problems, the corresponding action potential is nothing but the minimal 
action between two configurations for the Maupertuis action functional. This 
last can be defined as the energy functional associated to the Jacobi metric in 
the zero energy level. 

1.2 On the weak KAM theory. 

In order to give applications, we will show that theorem [5] enables us to prove 
a weak KAM theorem in the spirit of [lOj . HU. The novelty in this viewpoint, 
is that we regard the action of the Lax-Oleinik semigroup on a space of Holder 
functions. 
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Let us remember that a function u : —>■ R is said dominated by L, if 

it satisfies the condition u{x) — u{y) < (j){x,y) for all x,y € . Since the 

action potential is symmetric, theorem [5] implies that dominated functions are 
Holder continuous. On the other hand, it is not difficult to prove that they are 
locally Lipschitz in the open subset of total measure O C E^ ^ see proposition 
[7] below. Therefore, dominated functions are differentiable almost everywhere. 
We shall discuss this in more detail below. Another way to define the set 
of dominated functions, is using the Lax-Oleinik semigroup: given a function 
u : E^ —>• [— 00 , + 00 ) and t > 0 we define T^u : E^ [— 00 , + 00 ) by 

Tt~u{x) = inf { u{y) + (j){x, y,t) \ y & E^ ) . 

Then, a continuous function u is dominated if and only if u < for alH > 0. 
Notice that the set of dominated functions is convex and stable under the Lax- 
Oleinik semigroup. Setting = u for any function u, we will prove that 

{Tn t>o is a continuous semigroup on the set of dominated functions equipped 
with the topology of uniform convergence on compact subsets. 

Another set which is stable by the Lax-Oleinik semigroup is the set of 
functions which are invariant by symmetries. If we observe that the group 
of isometries of E, acts naturally on E^ by symmetries of the potential func¬ 
tion, then an obvious question is the existence of invariant fixed points of the 
semigroup. More precisely, we will say that a function u : E^ —> K is invariant 
if u(ri,..., rjv) = u(Ari + r,..., Arjv + r) for all x = (ri,..., tn) G E^, r G E 
and A G 0{E). 

Theorem 3 (invariant weak KAM). There exists an invariant and dominated 
function u : E^ —>■ M such that u = Tffu for all t > 0. 

In section[21 we prove the weak KAM theorem, and we study the relationship 
with the Hamilton-Jacobi equation. More precisely, we show that weak KAM 
solutions are global viscosity solutions in 17. 

An important difference with the compact case is that here the Aubry set 
is empty. In particular the technique used in m to prove the invariance of all 
solutions is not available. Moreover, we will exhibit non invariant solutions for 
the Kepler problem in the plane, which is the subject of the last section. 

2 Holder regularity of the action potential 

This section is devoted to the study of the action potential, and to give the 
proofs for theorems [1] and [7] 

2.1 Proof of theorem [H 

Given r G E and i? > 0, we say that a configuration x = (ri,..., r^) G E^ is 
contained in the ball B{r, R) when we have || — r || £; < i? for alH = I,..., iV. 

Suppose now that we have two configurations x and y such that for some r G E 
and some i? > 0, both x and y are contained in B{r,R). If we tried to bound 
4>ix, y, T) with the action of a linear path, then two problems arise. The first one 
is that the linear path can present collisions in which case the action is infinite. 
The second one is that, even if the linear path avoid collisions, the distance 
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between two given bodies can be arbitrary small for both configurations, hence 
the action can be arbitrary large. Both problems are solved in the following way: 
fix an intermediate configuration p with sufficiently large mutual distances, and 
take the linear path from x to p defined on [0,T/2] followed by the linear path 
fromp to y defined on [r/2, T]. This path has no more than 2N{N—1) collisions, 
and we can determine the values of t S [0, T] in which these collisions happen. 
Thus, reparametrizing the path in such a way that in the new times of collisions 
the action integral converges, we obtain the following proposition, from which 
we can easily deduce theorem [T] 

Proposition 4. Given two configurations x, y G contained in a ball B{r, R), 
r G E, R > 0, and given T > 0, there is a curve 7 G C{x,y,T), such that ^{t) 
is contained in B{r,(!>NR) for all t G [0,T], 

^ [ \'f{t)\'^dt<aT~^R^, and f UK,{'^{f)) dt < ljTR~'^'^, 

2 Jo Jo 

where a and ft are positive constants that only depend on the number of bodies, 
the total mass and the degree of homogeneity of the potential function. In fact 
we can take 

a = Qm\^MN'^ and /3 = 2 . 

1 — K 1 — K 


Proof. We first observe that it suffices to give the proof for a fixed value of 
T > 0: for ^ > 0, we can define a : [0, S'] ^ E^ as a(s) = jlsT/S), and we 
have 

nS pT 

/ I fT(s) I 2 ds = T^S-^ / I j{sT/S) I 2 ds = S-^T / I 7(t) | ^ dt < 2a R-^R^ , 
Jo Jo Jo 


S S T 

[ U^{a{s))ds= j U^{-i{sT/S))ds = ST-^ j U,,{jit)) dt < p SR' 


Jo Jo Jo 

We will then give the proof for T = 2. Take v € E such that || = 6R, and 

define p = (pi,... ,pn) G E^ by 


Pi=r + {i-l)v, i = l,...,N. 


Therefore, the configuration p is clearly contained in B{r, 6{N — 1)R). Notice 
also that the mutual distances pij = || Pi — Pj || _e of p are greater than 6R and 
smaller than 6{N — 1)R. 

Let now x = (ri,..., rw) be a configuration such that || — r || e < R 

for all J = 1,...,A^. We consider the curve : [0,1] E^, defined by 

Zxit) = x + 'if>x{t){p — x), where ifx '■ [0, 1 ] —[0, 1 ] is an increasing function, with 
V'a;(0) = 0 &od ifxO) = 1, to be determined. Our aim is to choose the function 
ifx conveniently, in order to obtain a bound of A{zx) which does not depend on 


Recall that if u and v are two vectors in a Euclidean space, and v Q, then 
we have, for all real number A, 


iiu+Aur 



) 2 


Ikl 


< u, u 
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As a consequence, we get 


I u + Au II > II u I 


A + 


< u,v > 


and the minimum of || m + Az; || is reached for 

^ < u,v > 

We will use the notation Uij = ri — Vj and vtj = {pt — pj) — (vi — rj) for 
i < j. Thus, the mutual distances of the configuration z^it) can be written 
dij(t) = II Uij + ipxit)vij II . Therefore, taking A = ipx{t), u = Uij and v = Vij in 
the above considerations, we deduce that each mutual distance dij(t) verifies 

dij(t) > II Vij II E I '4’x{t) - tij I > 4i? I tpxit) - tij I , 


where 


tij — 


<C Uij t Vij ^ ^ 


II ^..11 


2 

E 


Since || Uij || _e < 2i? and || Vij || e > 4i? for all i < j, we obtain that | tij | < 1/2 
for all i < j. Therefore, since then number of pairs (*,j) with 1 < i < j < N 
is bounded by iV^/2, assuming lemma [5] below, we can choose the function ijjx 
can be chosen in such a way that, on one side. 


'ipx{t)'^dt < 5iV 


21 + ^ 
1 — K 


and on the other side, for each i < j there is a real number Sij for which 


\i^xit)-Uj\ > N ^11 - Sij I 


for all t G [0,1]. Let us estimate the action A{zx) for this function tpx- We have 
Zx{t) = tpx{t) {p — x), and || Pi — ri\\ e < SNR, for alH = 1,..., TV. Hence, by 
the previous estimates we deduce that 


1 ^ 

/ \zx{t)\‘^dt = ■z'^mi\\pi-ri\\% j ipx{tfdt 

Jo ^ 1 ^0 


1 — K 


and that 


[ U^{zx{t))dt = V] /" rriimj dijit) dt 

Jo Jo 

Jo 
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Using that for r S (0,1), and for any s G R, we have 


\t- s\ 


■ dt = 


■ du < 2 


nl j 


du = 2 (1 — r) 


-1 


we deduce that 


[ U^{zS))dt < 
do 1 - « 

To finish the proof, let j/ = (si,..., s^v) be a second configuration contained 
in B{r,R), and define 7 € C{x,y,2) as follows: 7 (t) = Zx{t) if t < 1, and 
”f{t) = Zy(2 — t) if t > 1. We conclude that 


^( 7 ) = A{zx) + A{zy) < 320 


1 + K 


MN‘^ 


This also proves the proposition for T = 2, with 


4 iJl^ 7 VU«:+ 2 ) 7 Vf 2 j^-2k ^ 

1 — K 


a = 6Aol^MN'^ and /3 = 2 


1 — 


1 — 


It is important to note that for the Newtonian case (k = 1/2), the dependence 
of both constants in the number of bodies is in N'^. □ 

We have used the following lemma. 

Lemma 5. Given k € (0,1) and real numbers ai < ... < am, there are real 
numbers bi < ... < bm and an increasing absolutely continuous homeomorphism 
F 0 / [0,1] such that 


1 . 


\F{t)-ai \ > -^\t-bi\^/^+ 


2m 

for all t G [0,1] and each i = 1,... ,m, and 


2 . 


f F'{t)‘^dt < {4 +2a){m + 1) 
Jo 


1 + K 


1 — K ’ 

where a = min { | oi | ,..., | Om \ }• 

Proof. Given c > 0, let (/c : R \ { 0 } —5- M defined by gdx) = c | a; | 

Given b = {bi,. ■ ■ ,bm) € such that &i < ... < bm, we also define the 

function /f,_c : R \ { &i, ■ • ■, R by 

fb,cft) = niax{5c(i - &i), ■ • -jgdt -&„)}• 

We will define the required function U as a primitive of a function fb^c for a 
good choice of b and c. More precisely, we define F : R —>■ R. by 


F{t) = [ /b,c(s)ds- 

Jo 


The map F is an increasing homeomorphism of 
absolutely continuous. Moreover, we have 


As any primitive, F is 
\F{t)-F{bd \ >c(l + b:)|t-&,|i/i+" 
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Figure 1: Graph of fb,c for b = {bi,b 2 ,bz)- 


for all f G M and each i = 1,... ,m. Therefore, we must choose c > 0 and 
6 = ( 6 i,..., G K™ such that -F(l) = 1 and F{bi) = Ui for alH = 1,..., m. 

If we fix c > 0, then we have a unique possible choice for b. To see this, first 
observe that the m — 1 distances between the consecutive values of determine 
the m — 1 distances between the consecutive values of bi. If we set Ai = Oi+i — at 
and Bi = 6^+1 — bi then we must have 




/•bi+i /■Si/2 

/ fb,c{s)ds = 2c = + 

Jbi ’ JO 


hence Bi = 2 '^[Ai/c{l + From the condition F{bi) = oi we deduce 

that 

l>b\ r — bi 

ai= fbAs)ds = -j fb,c{s + h) ds . 

Jo Jo 

Therefore bi must be the unique solution of the equation ^ fb',c{s) ds = —ai 
where b' = {0, Bi, Bi + B 2 ,. ■., Bi + ... + Bm-i) = (0, 62 - &i, • ■ •, - &i)- 

Moreover, we have showed that there is a continuous vector b{c) G such that 
fo‘ fb(c),c(s) ds = Ui for all / = I, ..., m. Therefore, it is clear that 


6{c) 



also depends continuously on c. We claim that there is c G [l/2m(I + k), 2 + a] 
for which 6{c) = 1. We have 


m „i 

die) = /b(c),c(s) ds < V / gds-hic))ds . 

Jo Jo 


Since 


cl — bi (c 


gc{s — bi{c)) ds = / c\u\ < 2c / u = 2 c (1 +/^), 


/o 


' —bi{c 








Figure 2: Graph of F{t) = f* fb,c{s) ds. 


we deduce that (5(c) < 1 when c < l/2m(l + k). In order to prove the claim, 
it suffices to show that (5(c) > 1 when c > 2 + a. Since a = | Oj | for some 
j G { 1 ,..., m }, we have 


> 


> c 


bjic) 

ds I 

fbj(c) 

/ c\s-b,{c)\-^/^+^ 

Jo 
f\bAc)\ 

/ s-"/^+"ds = c(l + K)|6j(c)|i/^+", 

Jo 


' ds I 


which implies | hj{c) \ < [a/c(l + On the other hand we have 


(5(c) > mm{ 5 c(s - &j(c)) I s € [0,1] } 

> c(l + | 6 ,(c)|)-'‘/l+^ 

Thus, it suffices to prove that | hj{c) \ < _ i when c > 2 + a. By 

the previous estimation of | hj{c) \ , we only have to prove that (a /1 + < 

gi+K(ci+K/K _ condition is clearly satished if c > 2 and c> a. 

We take c G [1/2 to( 1 + «;), 2 + a] such that (5(c) = 1 and we define F : [0,1] —> 
[ 0 , 1 ] by 

Fit) = [ /b(c),c(s) ds . 

Jo 

In order to see that this function satisfy all the required conditions, it remains 
to estimate the norm of F'. If we observe that [0,1] \ { &i(c),..., &m(c) } has 
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at most m + 1 components Ij, and that on each one of these components we 
have 



we conclude that 


f F'{t)'^dt < (A + 2a){m + 1) 
Jo 


1 + At 

1 — At 


□ 


2.2 Minimal configurations. 

The following observations show that theorem [T] is optimal in the sense that the 
bound is reached by some configurations. We shall first recall the notions of 
centra] and minimal configurations, as well as some properties (see for instance 
Wintner [16], where the Newtonian case is discussed). 

We say that a configuration x £ is minimal, if it is a minimum of 
the potential function restricted to the sphere { y £ \ I{y) = I{x)}. 

In particular, minimal configurations are central configurations, that is, critical 
points of Uk. = I‘^ Uk, or in other words, configurations x £ which are critical 
points of Un restricted to { y £ E^ \ I{y) = I{x)}. Central configurations are 
also characterized as configurations which admit homothetic motions. In other 
words, a configuration xq £ E^ is central, if and only if Uk,{xo) < +00 and 
x{t) = r{t)xo is a solution of the N-body problem for some positive real function 
r(t). 

Take xq £ E^ a central configuration. If we look for an homothetic motion 
through Xq , then we must solve a one dimensional differential equation satisfied 
by r{t). A particular solution, that we shall call parabolic, is given by x{t) = 
Xf) for some value of c > 0. A simple computation shows that the action 
of this solution is 


Mx |[0.T]) 


2(1 + k)^ Jq 
2 


pT pT 

/ t-'^'^/^+^dt + c-2«C/«(xo) / 

Jo Jo 

Ji(1-k)/(1+k) ^ 


-2ktt / , 1 + «: 

2(1-ac2) 


If we set Rt = || x{T) || = || cxq ||, then we can write 

A(x|[o,T]) = aoT-iA?, + PoTRt'^\ 


for a good choice of constants ao and (3q. 

On the other hand, we will prove that if xg is a minimal configuration, then 
the above solution x(t) is globally minimizing. In other words, we have 


4'{0,x{T),T) = A{x |[o,T]), 

for all T > 0, therefore the bound for (l){x,y,T) given by theorem [1] cannot be 
improved modulo the choice of the constants. We will assume for simplicity 
that Xq is also normal, meaning that I{xo) = 1. 
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In order to prove that x is globally minimizing, we will first study the homo¬ 
geneous one center problem in dimension one which is satisfied by the function 
r{t) = More precisely, we have that r{t) must be an extremal for the 

Lagrangian system in R’*' = [0, -foo) defined by 


Lo{r,v) 


2 


where Uq = Uk{xo). 

It can be proved easily using the lower semi-coninuity of the Lagrangian 
action of Lq that, given 0 < ri < r 2 , there is at least one absolute minimizer 
between ri and r 2 . This means that there is a curve 7 (ri,r 2 ) : [OjT'] R+, 
for some positive time T > 0, such that 7 (r’i,r 2 )( 0 ) = ri, 7 (ri,r 2 )(T) = r 2 , 
and such that 7 (ri,r 2 ) minimizes the Lagrangian action in the set of all ab¬ 
solutely continuous curves a : [a,b] 1R+ with a{a) = ri, a{b) = r 2 , and 

a < b. Moreover, the fact that this absolute minimization property holds 
(i.e. with fixed extremities but in free time), implies that the energy of the 
extremal 7 (ri,r 2 ) must be critical (zero). Therefore 7 (ri,r 2 ) satisfies the dif¬ 
ferential equation 7 ^ = 2 [/o 7 “^^^ and we conclude the uniqueness of such ab¬ 
solute minimizer. By integration we get that for any r > 0, the absolute mini¬ 
mizer 7 ( 0 , r) : [0,r] —)■ R+ is defined for T = (I -|- and that 

7 ( 0 , r){t) = where c = -|- 

We prove now that if xq is a minimal configuration, then the parabolic 
motion x(t) = defined above is globally minimizing. Fix T > 0, and 

take any other curve 7 G C(0, x(T), T). We must to prove that ^( 7 ) > A{x |[o_t] 
). In fact we will prove that if S' = sup {t G [0, T] | ^{t) = 0 } then 

A{l |[S,T]) > A{x |[0,T]) • 

Setting 7 i = 7 |[5 y] we have 71 (t) 7 ^ 0 for all t G (S, T]. Thus we can set 
7 i(t) = r{t)s{t), where r{t) = | 7 i(<) | = /( 7 i(t))^/^ for t > S. Obviously we 
have that | s{t) \ = I for all t > S, and the action of 71 can be written 

^(71) = S / ^ \ f ^dt + f r{t)-^^U^{s{t)) dt. 

J S J S J s 

Since xq is minimal, we have U(s{t)) > Uq = Uk{xq) for all t > S. Moreover, 
we have ^ 

^( 71 ) > A{rxo) = y Q f{tfdt + C/or(t)“^'"^ dt. 

Note that the last integral is nothing but the Lagrangian action of the curve 
r{t) : [S, T] R’*' for the one dimensional homogeneous one center problem, 
and as we have see the absolute minima of the Lagrangian action is reached by 
the zero energy solution r{t) = But in that case, the last integral is 

also de Lagrangian action of the parabolic motion x{t) restricted to the interval 
[0,r] for the N-body problem. We deduce that 

^( 7 ) > ^( 71 ) > y f{tfdt -b [/or(t)-2«^ dt = A{x |[o,T]) 
hence we conclude that the solution x{t) is globally minimizing. 
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2.3 Properties of the action potential and proof of theo¬ 
rem [21, 

We start this section by showing that the action potential is a distance function 
on . 

Proposition 6 . For all x,y G we have (j>{x,y) = 0 if and only if x = y. 

Proof. Let x G E^ be a configuration, and choose a path a : [0,1] —?> E^ which 
satisfies cr(0) = x and A{a) < + 00 . Then define for 0 < T <2 the curve 
7 t G C{x,x,T) by ■jTft) = o-(t) if t < T/2, and 7 T(t) = cr{T — t) if t > T/2. 
It is not difficult to see that -g 0 as T ^0, from which it follows that 

(p{x, x) = 0 for all x G E^. 

To see that the condition is necessary, take any two configurations x = 
(ri,... ,rAr) and y = {si,.. .,sn) in , and a path 7 = ( 7 ^,.. .,jn) G C{x,y). 
If d = II y — a; II , and 7 is defined on [0,r], then it must exist Tq £ [Oj^] such 
that II 7 (To) — a: II = d and || 7 (t) — x || < d for alH e [0, Tq] . Moreover, we must 
have d = || 'fiiTo) — ri\\ e for some f G { I,..., iV }. If Tq > I, we can write 

A( 7 ) > d (7 |[o,To]) > [ dt > C > 0 , 

Jo 

where C = min { Uk,{z) \ || z — x || < d }. If Tq < 1 we have 

. / ^ w I X m-i II . X X II 9 , TO d^ 

d(7) > d (7 |[o,ro]) II \\Edt> , 

where to = min { toi, ..., ijin }■ The last inequality follows from the fact that 
7 i is absolutely continuous and the Cauchy-Schwartz inequality. Therefore, we 
conclude that if (j){x, y) = 0 , then d = 0 and x = y. □ 

In the sequel we will denote <5(z) the minimal distance between the bodies of 
the configuration z. More precisely, S : E^ —>• R+ will be the function defined 
by 5{z) = min { || Zi — Zj || e | i < j }, where z = (zi,..., zn)- Thus the set 
of configurations without collisions is nothing but = { z G E^ \ i5(z) > 0 }. 
From the following proposition we can easily deduce that the action potential 
is a locally Lipschitz function in 17 x 

Proposition 7. There are continuous functions k{z) > 0 and e(z) > 0 m 17 
such that, if x,z G E^ satisfy || x || < e(z), then (j)(z, z + x) < fc(z)|| x || . 

Proof. We give the proof for e(z) = d(z)/4. Since z is without collisions, we have 
e(z) > 0. For T > 0 we define the curve 7 : [0, T] — 5 > E^, by 7 ( 7 ) = z + {t/T)x. 
If z = (zi,...,zjv) and x = (ri,..., xat), then 7 ( 7 ) = ( 71 ( 7 ),..., 7 Ar( 7 )), where 
7 i( 7 ) = Zi + {t/T)ri. Hence, for i < j and 7 G [0,T] we can write 

7ij W = II %it) - Ijit) II E > II - Zj II E - (7/r)|| ri - rj || e > S{z)/2 , 

and 

fd„(7(7)) = ^ruiTOj 7ij(7)“^'‘ < [d(z)/ 2 ]“^'^ . 

i<j 
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Therefore, using that | x | ^ = I{x) < MN\\ x || we deduce that 
^( 7 ) = ^ \x/T\^dt + U{-i{t))dt 

< MiV||x|| V2T + M^iV^ [5(z)/2]"^” T. 

If X = 0 there is nothing to prove, since we already know that (j){z, z) = 0. If 
X ^ 0, we can take T = || x || , and the above estimation gives ^( 7 ) < k{z)\\ x || 
for k{z) = MN/2 + {5{z)/2)~^^. □ 

We introduce now a notion of cluster partition of a subset A C E adapted 
to our purposes. Given A > I, we will say that the set { ri,..., r/f } C i? defines 
a A-cluster partition of size i? > 0 of 7l, if the following two conditions are 
satisfied: 

1. II Ti — rj II E > 2Ai? for all 1 < * < j < K, 

2. A is contained in the union B{ri,R). 

It is clear that if A is finite and R is small enough, then A defines itself a cluster 
partition of size R of A. It is also clear that if A is bounded, then any r G A 
defines a trivial cluster partition of size R for any R > diam(A). 

We will need the following lemma. 

Lemma 8 . Given A > 1, A = { ri,..., rjv } C i? and e > 0, there is a subset 
A' C A, and R{e) > 0 such that: (i) e < i?(e) < (2A)^e, (ii) A' defines a 
cluster partition of size i?(e) of A. 

Proof. We reason recursively. We begin setting A[ = A. If A[ does not define a 
A-cluster partition of size e, then there are r,sG A'^ such that || x — s || £; < 2Ae . 
If that is the case, we define ^2 = \ { s }. Then we reason as before: if A '2 

does not defines a A-cluster partition of size 2A e then we have r, s G A^ such 
that II r — s II < ( 2 A)^ e , and we set A 3 = ^2 \ {■^ }• It is clear that the process 
finish at the most in N steps. □ 

Using the existence of cluster partitions we will prove the following proposi¬ 
tion, from which theorem [5] can be deduced as a simple corollary. 

Proposition 9. There are positive constants oi, /3i >0 such that, for all 
x,y G and for all T > 0 we have 

(^{x,y,T)<aiT-h^ Te”"", 

whenever e > || x — y || . The constants oti and (3i only depend on the degree of 
homogeneity of the potential (—2 k), the number of bodies N, and their masses. 

Proof. Fix a configuration x = (xi,... ,rjv) G E^, and denote by A^ the set 
{ xi,..., r^r } C E. Let y = (si,..., sw) G E^ be any other configuration. If 
we apply lemma |S] to A^ with e > || y — x || and A = 247V, we conclude that 
there are ,..., x^j^ G A^, and R(e) > 0 with the following properties. 

1. e < R(e) < (487V)^e, 
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2. for all 1 < j < fc < if, we have || || e > 48A^i?(e), and 

3. A^UAy is contained in the disjoint union Uj=i where Bj = B( ri. , 2 i?(e)). 

Therefore, both configurations x and y are decomposed in K clusters, each 
one contained in a ball Bj. More precisely, we have a partition { 1,..., iV } = 

Ii U ... U Ik such that i S Ij if and only if both and Si are in Bj . Denote 
by Nj = card(/j) the number of bodies in cluster j, and by Mj the total mass 
of this cluster, that is Mj = Thus we have N = Ni + ... + Nk and 

M = Ml + ... + Mk- 

We consider now the Nj -body problem composed by the bodies in the ball 
Bj . Given T > 0, we apply proposition |4] in each ball Bj , j = 1,... ,if, with 
initial and final condition conformed by the Nj bodies of x and y contained in 
Bj . Therefore we obtain, a path 7 = ( 71 ,..., jk) G C{x, y, T) such that for all 
j = 1 ,..., if we have, 

1. If i e Ij , then ^i{t) G B{ri., 12N R{e)) for all t G [0,T], 

2 . 

= \ \\Edt< 10® ^ M, Nf R{ef/ T , and 

do , 6 /^ ^ 

3. 

Wj= ^ TOiTOfcll 7iW - 7fcW II 
do i.keij 

< 2^^^ N^'^+^Mf 12-^'^ R(e)-'^^T . 

I — K 

Notice that the action of the curve 7 = ( 71 ,... ,Jn) G C{x,y,T) is 

K K 

i=i i=i 

where Wq is the integral of the terms of the potential function 11^ corresponding 
to pairs of bodies in different clusters. More precisely, 

Wo= f ^ ^ m.TOfc ||7i(t)-7fc(t) ||s^'"dt. 

do i<j<i<K ieij,keii 

Since the balls B{ri.,24:NR{e)) are disjoint, we deduce that 
Wo < N^M^ (24iV)“^'" i?(e)"2'= T . 

Using that using that R{e) < {2AN)^ e, we can write 
A(7) < ai T-^e^ +/3i re- 2 « , 

for some positive constants ai and /3i only depending on N, M and k. □ 
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Corollary 10. There is a positive constant p, > 0 such that for all x G 

Proof. It suffices to take e = in proposition [S] □ 

Proof of theorem\^ By the previous corollary we know that 4>(x,x) — 0 for all 
X S . If x, 2 / € E^ are two different configurations, then proposition [5] says 
that 

<t>{x,y,T)<ai , 

for all T > 0 and for any e > || a; — y || . Since the right hand of this inequality 
is a continous function of e > 0 we also have 


(f{x,y,T) < aiT '^\\x-y\\'^ + /3iT\\x-y 


and the proof is achieved taking T = || a: — ?/1| 


□ 


2.4 Homogeneity of the action potential. 

There is a property of homogeneity of the action potential due to the homo¬ 
geneity of the potential function Uk- We did not use this property in the above 
proofs, but we think that it is useful to complete the picture of the action po¬ 
tential. The proof can be done reparametrizing conveniently homothetic paths 
of a given path. 

Proposition 11. If X> 0, then cj){Xx,Xy) = X^~^(j)[x,y) for all x,y € E^. 


3 Weak KAM theory 

The relationship between global solutions of the Hamilton-Jacobi equation and 
globally minimizing solutions of the corresponding Lagrangian flow is well known. 
Let us recall that the Hamiltonian, defined on T*E^ = E^ x [E*)^ is the func¬ 
tion 

Hix,p) = i |p|2 - U^{x ), 

where | p \ denotes the dual norm of p G {E*)^ with respect to the norm on E^ 
induced by the mass scalar product. More precisely, if we identify the space E 
with its dual E* using the scalar product < , >e, and p = (pi,...,pjv) G (E*)^, 
then 

N 

\p\^ = \\p^ III. 

A closely related function is the total energy, defined on TE^ as £ = H o 
C, where C : TE^ —>• T*E^ is the Legendre transform C{x; vi,...,vn) = 
(x] pi,. .. ,pn), Pi = TTiiVi. It is easy to see that f is a first integral of the 
motion. 

We will prove the existence of critical global (weak) solutions for the Hamilton- 
Jacobi equation H{x,dxu) = c. The critical value of this Hamiltonian can be 
defined as the infimum of the values of c G ffi. such that the Hamilton-Jacobi 
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equation admits global subsolutions. Since inf^w Uk,{x) = 0, and constants 
functions are global subsolutions for c = 0, it follows that the critical value is 
c = 0. Therefore, we are interested in global solutions of 

\d,u\^ =2U^(x). (HJ) 

We will obtain global solutions as fixed points of a continuous semigroup acting 
on the set of weak subsolutions, namely the Lax-Oleinik semigroup. There are 
no new ideas in the method that we apply here. In fact, we will follow the 
scheme introduced by Fathi in m. with some adjustments to our setting. As 
we have said in the introduction, the difference is that we consider a space of 
Holder functions on which the semigroup acts, and theorem [2] will assure that 
the method works with this space. 

3.1 The Lax-Oleinik semigroup. 

Given a continuous function u : —5> R and t > 0, we define T^u : —>■ 

[—oo, +oo) by 

It” u{x) = inf { u{y) + (t){x, y,t) \ y & E^ ) . 

We also define T^u = u for all function u. The semigroup property follows from 
the definition . In other words, for any function u we have that T^{T~u) = 
for all t, s > 0. We will restrict the semigroup to the set V. of dominated 
functions. More precisely, we define 

H = { u : E^ —>• R I u{x) — u{y) < (j){x, y) for all x, y € E^ } . 

Notice that u : E^ —)■ R is in H if and only if u < T^u for all t > 0. On 
the other hand, u < v implies that T^u < T^v for all t > 0. Therefore, the 
semigroup property implies that Tfu € H for all u € 7^. Also notice that 'H is 
convex, and nonempty since it contains all constant functions. 

In the sequel, the set Ti. will be endowed the compact open topology, that is 
to say, the topology generated by the sets 

Uk(u, e) = { X S 77 I I v{x) — u{x) \ < e for all x € TF } , 

with u € K, K C E^ compact, and e > 0. 

Proposition 12. The map T~ : 77 x [0, +oo) —>■ 77, (rt, 7) i— T^u is continuous. 
We will use the following lemma. 

Lemma 13. For all x, y G E^ and T > 0 we have <j){x,y,T) > (m/2T)|| x — 
y II w/iere m = min { TO 1 ,..., m^v }• 

Proof. Let r, s G E and a : [0, T] —^ if an absolutely continuous curve such that 
cr(0) = r and cr{T) = s. We observe that 
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hence 


\\r-s\\l<T [ \\a{t)\\ldt. 

Jo 

If a: = (ri,..., rjv) and y = (si,..., sjv) are two configurations, then we can 

choice i G { 1 ,..., } such that || — Si || _e = II a; — y || . Take now 7 = 

( 71 ,, 7 Af) G C{x, y, T). By the previous observation we have, 

^( 7 ) > i'mi/2) [ II 7 i(t) III > imilTT) || - 5^ || | > (m/ 2 r) \\x-y\\‘^ , 

Jo 

which proves the lemma since 0(x, y, T) = inf { ^( 7 ) | 7 G C(x, y, T) }. □ 

Proof of vrovosition[TR As a first step, we show that given R > 0 and t > 0, 

there is a constant k{R,t) > 0 such that 

Tt~u{x) = inf {u{y) + (j){x,y,t) \ || y - x || <k{R,t)} 

for all u G and all x G with || x || < R. To see this, fix i? > 0, t > 0, 

u € PL and x G such that || x || < R. Suppose that y G E^ is such that 
II y — X II > 1 and M(y) + 0(x, y, t) < u(x) + (/)(x, X, t). Then, by lemma [T3] and 
theorem [5] we have 

^ II y - a; II ^ < y II y - a; + <^(x, X, t). 

Therefore, using that || y — x || >1 and theorem [1] we deduce 

m II y — X II ^ < 2?7t II y — X II + 2aR^ + 2/3 t^R~^'^ , 
hence that || y — x || < ko{R, t) where 

ko{R,t) = rjt/m + /mf + 2aP?jm + . 

Setting k{R,t) = max { 1, fco(i?, t) }, it follows that u(y) + (j){x,y,t) > u{x) + 
(j){x, X, t) for all y G E^ such that || y — x || > k{R, t), and we conclude that 

Tf~u{x) = mi{u{y) + (j){x,y,t) \ y€E^) 

= inf {■u(y) + (/()(x, y, t) | || y - x || <k{R,t)} . 

Let now u,v G PL and t > 0. Let K C E^ be a compact subset, and R > 0 
such that II X II < R for all x G AT. If we set 

Kt= IJ { y G I II y - a: II < k{R, t) } , 

a:e K 

then for dX\ x ^ K we have 

Tf v(x) = inf { v(y) + <f(x, y,t) \ y G Kt} . 

On the other hand, since 

v{y) < u{y) + sup {I u{y) - u(y) I | y g a:* } 
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for all y S Kt, we deduce that T^v{x) < inf { u(j/) (f){x,y,t) \ y & Kt} + 

sup { I u(jj) — v{y) I \ y € Kt}. Thus we have proved that Tfv{x) — T^u{x) < 
sup { I u{y) — v{y) \ | y G itT* } for all a; G K. Hence we have that 

I v{x) - Tfu{x) I < sup { I u{y) - v{y) \ \ y & Kt} 

for all x G K. Since the subset Kt C is compact, this implies the continuity 
of the map for each < > 0. It remains to prove the continuity of with 
respect to t. Since {T^)t>o is a semigroup it suffices to prove the continuity at 
t = 0. Given u G "H, we have by corollary (TU] that 

0 < T^u{x) — u{x) < (j){x,x,t) < 

for all X G . Therefore, T^u converges uniformly to u when t —>■ 0. As 
we have said, using the semigroup property we can deduce that T^u converges 
uniformly to Tf~u when t ^ tg for alHo > 0. □ 

3.2 Proof of theorem [31. 

Proof. Let H be the quotient space of H by the subspace of constants functions. 
Thus, is homeomorphic to Ho = { u G H | m( 0) = 0 }. By theorem [5] we 
have that dominated functions are uniformly equicontinuous. It follows that 
Ho is compact by Ascoli’s theorem. Therefore, H is a compact, convex, and 
nonempty subset of C^{E^ the quotient of the vector space C^{E^ ,R) by 
the subspace of constant functions. Notice that R) is endowed with the 

quotient topology of the compact open topology on R). In particular, 

R) is a locally convex topological vector space. 

Since Tff [u + c) = Tffu + c for all c G R, it is clear that the semigroup 
T~ defines canonically a continuous semigroup Tff : H —H. If we apply 
the Schauder-Tykhonov theorem, see [8] pages 414-415, we conclude that Tf 
has a fixed point in H. That is to say, there is a function u G H such that 
Tffu = u -\- c{t) for some function c : [0, +oo) —^ R. The semigroup property 
and the continuity of T~ imply that c{t) = c(l)t. Since u G PL, we have that 
u < Tffu for all t > 0, hence we must have c(l) > 0. We will prove that 
c(l) = 0. Notice that Tffu = u + c{l)t implies u{x) — u{y) < (j){x,y,t) — c{l)t 
for all x,y G E^. Hence, by theorem [T] we have that 

u{x) - u{y) <a^+ - c(l)^ t 

whenever x and y are contained in a ball of E of radius R > 0. Since this 
must be true for R and t arbitrary large, we conclude that c(l) = 0. Therefore 
Tfu = u for all t > 0. 

It remains to prove that there are fixed points of T~ which are invariant 
by the group of symmetries. This can be done as in m as follows. We de¬ 
fine the Hint, as the set of functions in Pi which are invariant by symmetries. 
Thus PLinv is also convex, closed and nonempty since constant functions are 
invariant. Moreover, PLinv is stable by the Lax-Oleinik semigroup. Therefore, 
the quotient of this set by the subspace of constants functions is also compact, 
convex, nonempty and stable by the induced semigroup T~. With the same 
arguments as above we obtain an invariant fixed point. □ 


18 



3.3 Viscosity solutions and subsolntions. 

It is well known that the notion of dominated function is related to a notion of 
subsolution of the Hamilton-Jacobi equation, namely the notion of viscosity sub¬ 
solution. On the other hand, viscosity solutions (see below) can be detected as 
fixed points, modulo constants, of the Lax-Oleinik semigroup. An introduction 
to the subject of viscosity solutions can be found for instance in the books [T], 
[2] or [^. However, our setting presents some technical differences, essentially 
due to the fact that the potential function is infinite in the set of conhgurations 
with collisions. The following is a little adaptation of some results in section 5 
of [H]. 

Recall that u : —>• K. is a viscosity subsolution at Xq € of (HJ), if for 

each function ?/; : E^ -y R such that Xg is a maximum ot u — ijj we have 
I I ^ < 2 {/^(xo). Given V C E^, we say that u is a viscosity subsolution 
in V if it is viscosity subsolution at each x € V. We remark that any function 
is trivially a viscosity subsolution in where H C E^ denotes the set of 
configurations without collisions. 

Analogously, a function u : E^ —>• K is said to be a viscosity supersolution 
at xq e E^ of (HJ), if for each function if} : E^ R such that xg is 
a minimum oi u — ip we have | dxgip | ^ > 2 UK,{xg). If xg G then u is a 
viscosity supersolution at xg if and only if there are no functions ip such 
that Xg is a minimum ot u —ip. As for subsolutions, given V C E^, we say that 
M is a viscosity supersolution in V if it is viscosity supersolution at each x G V. 

We say that a continuous function u : E^ —>■ R is a viscosity solution of 
(HJ) in H C E^ if it is both a subsolution and a supersolution in V. It is not 
difficult to see that a viscosity solution u satisfies (HJ) at each point x G V 
where the derivative dxU exists. We will prove the following. 

Proposition 14. (1) Any u G 71 is almost everywhere differentiable and a 
viscosity subsolution of Hamilton-Jacobi in E^. 

(2) If u G H is a fixed point of the Lax-Oleinik semigroup, then u is a viscosity 
solution of Hamilton-Jacobi. 

Proof. The fact that dominated functions are differentiable almost everywhere 
follows from proposition [71 the fact that collisions are contained in a finite 
number of affine subspaces and the Rademacher’s theorem. In order to prove 
that they are viscosity subsolutions, take u G % and ip : E^ —>■ R of class 
such that u — Ip admits a maximum at some xg G E^. Let v G E^. For all 
t > 0 we have 

1 p{xg) — 1 p{xg — tv) < u(xg) — u(xg — tv) — 2 J J ^k(^0 + Sv) ds . 

Dividing by t and taking the limit for t —>■ 0 we obtain 

dxo'P{v) < PU^,{xg) . 

If we define pi,... ,pn G E by the condition dxgipiv) = < Pit^i >e for 

all V = (wi, ..., vm) G E^, then we can write 

N 

\dxo'<P\^ = Y1 II P* II i = dxoi’iw) 


19 


where w = (wi ,..., wn) and rriiWi = pi for all z = 1,..., iV. Therefore, using 
the last inequality with v = w we obtain | dxQtp I ^ < 2 Uh{xq). 

It remains to prove (2). Suppose that u € H is such that T^u = u for all 
< > 0. Let ijj : —>■ R be a function such that u — -tp has a minimum at 

some xq € fl. 

With the same arguments as in the proof of proposition IT^ we deduce that 
there is a constant k > 0 such that 

Tf u{xo) = inf { u{y) + (j){xo, y,l) \ || y - a;o || < k} . 

Therefore, using theorem [1] and the lower semi-continuity of the Lagrangian 
action we can choose yo € such that || ?/o “ II < ^ and a curve 7 G 
C(a;oi?/ 0 )l) such that 

u{xo) = T~u{xo) = u{yo) + ^ f \'p{t)\‘^dt+ f U^{'^{t))dt. 

^ Jo Jo 

In particular, since m is a dominated function we must have 

u{xo) - u{j{t)) = ^ [ \j{s)\‘^ds+ f U^{-f{s))ds 
Jo Jo 

for all t G [0,1], which says that 7 is a calibrated curve for u. Using that a;o is 
a minimum for u — ip we conclude that 

ip{xo) - ip{l(t)) > u(xo) - u{-f{t)) = \ j^\ 7(s) \^ds + U^i-fis)) ds . 

At this point we can prove that xq is a configuration without collisions. In 
other words, u is a viscosity supersolution at each collision configuration xq 
because there are no test functions ip such that u — ip has a minimum in xq. 
To see this we proceed as follows: 

Let AT > 0 be a Lipschitz constant for the restriction of ip to some neigh¬ 
borhood of xo- Using the fact that 17^ > 0 in the last inequality we deduce 
that ^ 

a: |7(t) - 7 ( 0 ) \ l7(s) l^'^s 

for any f > 0 small enough. Applying the Cauchy-Schwarz inequality we also 
have ^ 

| 7 (s)|^ds > | 7 (s)|ds^ >|7(i)-7(0)r 

hence 

|7W-7(0)| <2Kt. 

Using again the Lipschitz constant for ip and the previous inequality, but ne¬ 
glecting this time the kinetic term, we obtain 

2tK‘^ > [ Ut^{'y{s))ds . 

Dividing by t and taking the limit for t ^ 0 we deduce that Uk,{xq) < 2K^, 
meaning that xq G D. Since xq is a configuration without collisions, there is 
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<5 > 0 such that 7 ([ 0 , 5 ]) C ri. Since 7 is minimizing, hence a solution of the 
Euler-Lagrange flow in [0, 5] we know that 7 is differentiable at t = 0. Dividing 
by t and taking the limit for t —>■ 0 in the inequality 

ipixo) - tpilit)) > ^ \j{s)\^ds + J^ C/„( 7 (s))ds 

we obtain 

dxo^l^iv) > 2 + 

where v = — 7 ( 0 ). On the other hand, always we have 2p{v) < |p | ^ + | f | ^ for 
p G (E*)^ and v G . Thus we conclude that | dxoi’ | ^ > 2 Uk{xo). We have 
proved that m is a viscosity supersolution at Xq. □ 

3.4 Lax-Oleinik and weak KAM solutions. 

Following the analogy with the weak KAM theory for Tonelli Lagrangians on 
compact manifolds, we show that the fixed points of the Lax-Oleinik semigroup 
are the weak KAM solutions defined by Fathi in m- More precisely, we show 
that the fixed points of Lax-Oleinik semigroup are characterized by the following 
property: given any configuration x G E^, always we have a calibrated curve 
jx ■ (— 00 , 0 ] —>■ E^ such that 7 a; ( 0 ) = x. 

Recall that if u : E^ ^ R is a dominated function, then a curve 7 : / ^ E^ 
is said to be calibrated when satisfies 

m( 7(6))-u(7(a)) = A(7|[„_6]) 

for all compact interval [a, b] C I. In particular, the calibrated curves of a 
dominated function are free time minimizers, meaning that 

4(7l[a.h]) = 'i^’(7(a),7(&)) 
for all compact interval [a, b] C /. 

Therefore, the fixed points of the Lax-Oleinik semigroup can be characterized 
in terms of calibrated curves as follows (recall that our Lagrangian is symmetric). 

Proposition 15. Let u G TL be a dominated function. Then u = Tffu for all 
t > 0 if and only if, for each x G E^ there is a curve 7^; : [ 0 , -I-00) —>■ E^ with 
7a;(0) = X and such that u(x) = u{'jx{t)) + A( 7 |[o_t]) for all t > 0. 

Proof. Suppose first that the condition is satisfied. Take x G E^ and the 
corresponding calibrated curve ^x '■ [0,-l-oo) —>■ E^ with 7 a; (0) = x. Since 
u G % already we known that u < Tffu for all t > 0. On the other hand, if we 
fix t > 0 we have T^u{x) < u{pjx{t)) -I- A( 7 |[q jj) = u(x). Therefore u is a fixed 
point. 

Suppose now that u G TL is a fixed point of (J)~). Given a configuration 
X G E^ and t > 0 we have 

uix) = T~u{x) = inf { u{y) + (j)[x, y,t) \ y G E^ ] . 

Using lemma [T3] and theorem [5] (as in the proof of proposition [T^ we deduce 
that there is a constant k > 0 (depending on x and t) such that 

u{x) = Tfu{x) = inf { u{y) + (j){x, y,t)\yGE^and\\y-x\\ <k} . 
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Therefore, using theorem [T] and the lower semi-continuity of the Lagrangian 
action we can choose y{x,t) G such that || y{x,t) — a; || < k and a curve 
Jx.t G C{x,y{x,t),t) such that 

uix) = T['u{x) = u{y{x, t)) + 

For each positive integer n > 0 we define the curve 7 „ : [0,n] ^ as the 
curve 7a;, ra- Observe that if m > n then 7^ |[o,ri] minimizes the action in 
C(x,^m('n-),n). Now we apply theorem [T] and once again lemma [T51 and we 
deduce that for a fixed positive integer n > 0 , the sequence (A( 7 m |[o,r!,]))m>n 
is bounded. It is not difficult to see (using the Cauchy-Schwartz inequality) 
that an absolutely continuous curve 7 : / —>■ E^ with finite Lagrangian action 
must satisfies | 7 (t) — 7 ( 3 ) | < 2A('y) 11 — s | for all t,s € I. Then we can 
apply Ascoli’s theorem and deduce the existence of a convergent subsequence 
of ( 7 m |[o.ra])m>n- By a diagonal process we can extract an increasing sequence 
of indexes G N such that, for each positive integer n > 0, the sequence 
ilrrik \[o,n])mk>n converges uniformly, when k — ^ 00. Observe now that by 
construction, each curve qm^ |[o,n])mfc>n calibrates the function u. Therefore 
the semi-continuity of the action implies that the curve '■ [ 0 , -l-oo) — E^ 
defined by 7 x(t) = limfc_>oo 'Ymkit) is also calibrated. □ 

We remark that for the Newtonian potential (k = 1/2), Marchal’s theorem 
implies - except of course in the collinear case (dim A = 1 ) - that the calibrated 
curves of weak KAM solutions are true motions for t > 0 since they must be 
contained in 0, the set of configurations without collisions. The dynamics of 
the free time minimizers of the Newtonian N-body problem is described in [7]. 


4 The Kepler problem 

Unfortunately, the proof of the weak KAM theorem do not give any explicit 
solution. Nevertheless, we can give explicit solutions for the Hamilton-Jacobi 
equation of the Kepler problem. 

We will find first isometry invariant solutions when we have two bodies of 
unit mass in a line {N = 2, mi = m 2 = 1, k = 1) and a Newtonian potential 
(k = 1/2). An invariant solution u : —>■ R must satisfy u{x+z, y+z) = u{x, y) 

and u{x, y) = u{—x, —y) for all x^y, z G M. Therefore the solution must be of 
the form u{x, y) = f{\x — y\), and the Hamilton-Jacobi equation reads 

ul + ul = 2\x-y\-^ . 

Replacing u{x,y) by f(\x — y\) and solving the differential equation in / we 
conclude that the unique invariant global solutions (up to an additive constant) 
are the functions 

u±ix,y) = ±2\x . 

In fact, the positive solution is the unique invariant fixed point of the forward 
Lax-Oleinik semi-group 

T+uix) = inf { u{y) - (j){x, y, t) | y G } 

and therefore, the negative one is the unique invariant fixed point of the back¬ 
ward semigroup T['. Of course, since the Lagrangian is symmetric in speed, we 
have that u G 7^ is a backward solution if and only if —u is a forward solution. 
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Figure 3: The two solutions and their derivatives for the 1-dimensional Kepler 
problem. 


It is not difhcult to see that we have also solutions invariant under transla¬ 
tions, in particular the function 6_|_ given by b+{x,y) = U-{x,y) ior x > y and 
b+{x,y) = u+{x,y) tor x < y is also a weak KAM solution. 

For the planar Kepler problem it is convenient to reduce first the problem 
by fixing the center of mass at the origin, or equivalently, to look for translation 
invariant solutions. Since the configuration is then determined by the position 
of the first body a: € the problem reduces as usual to the center fix problem. 
If we denote x = {xi,X 2 ) the position of the body, then the Hamilton-Jacobi 
equation reads 

+ ul^{x) = 2\\x\\ . 



Figure 4: Calibrated curves of solutions for the planar Kepler problem. 
Doubtlessly, the simplest solution that we can give is the rotation invariant 
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solution 


u{xi,X2) = -{xf + . 

Its calibrated curves are all the parabolic homothetic motions, represented in 
the left side of figure 21 The half parabolas at the right side are the calibrated 
curves of a Buseman type solution, which is constant and not differentiable over 
the dashed line. A computation made by A. Venturelli shows that this last 
solution can be explicitly defined by the formula 

u{xi,X 2 ) = -((xf + x\Yl'^ + . 
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